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Abstract 

Inhomogeneous quantum cosmology is modeled as a dynamical system of discrete 
patches, whose interacting many-body equations can be mapped to a non-linear mini- 
superspace equation by methods analogous to Bose-Einstein condensation. Compli- 
cated gravitational dynamics can therefore be described by more-manageable equa- 
tions for finitely many degrees of freedom, for which powerful solution procedures are 
available, including effective equations. The specific form of non-linear and non-local 
equations suggests new questions for mathematical and computational investigations, 
and general properties of non-linear wave equations lead to several new options for 
physical effects and tests of the consistency of loop quantum gravity. In particular, 
our quantum cosmological methods show how sizeable quantum corrections in a low- 
curvature universe can arise from tiny local contributions adding up coherently in 
large regions. 



1 Introduction 

One of the main problems in deriving a reliable Planck-regime scenario in canonical quan- 
tum cosmology is the question of how to include inhomogeneity. While homogeneous 
models can easily be quantized, inhomogeneous degrees of freedom severely complicate 
mathematical evaluations. Even the formulation of consistent evolution equations, sub- 
ject to the anomaly problem, remains incomplete: no consistent and covariant version of 
inhomogeneous modes valid at high density and including all relevant quantum effects is 
available at present. 

As possible solutions, two approaches have been developed so far, mainly with the 
methods of loop quantum cosmology [HE]- First, effective equations have been successful 
in addressing the anomaly problem [31 H] and in including all relevant quantum effects in 
sufficiently general form. (For details, see [5].) Potentially observable phenomena have 
been uncovered, showing physical consequences of discrete quantum geometry [U El [7J El [9j 
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[TUl [TT| [T2| IT3] and making the theory falsifiable [TJ1 [15] . At high density, the implications of 
quantum space-time can be dramatic: general properties of effective constrained systems 
show the presence of signature change, turning Lorentzian space-time into a quantum 
version of 4-dimensional Euclidean space [HI \TT\ [T8j EE]. Bounces as they had often been 
envisaged as non-singular versions of cosmology [20], and formulated in quite some detail 
[21] in homogeneous models of loop quantum cosmology, are then replaced by acausal 
pieces of 4-dimensional space devoid of deterministic evolution. Regarding specific field 
equations and details of the transition, however, present calculations remain incomplete 
because not all quantum effects could yet be implemented consistently at high density 
Moreover, although state properties can be derived by canonical effective equations, finding 
full quantum states is difficult in this setting. 

One of the alternatives is so-called hybrid quantization [22], in which one combines 
a loop-quantized homogeneous background model with Fock quantized inhomogeneous 
modes. Wave functions can then be solved for and evolved, at least numerically and with 
certain truncations [23] • However, using a Fock quantization for inhomogeneous modes, one 
does not directly deal with the discreteness of space-time. Moreover, the hybrid method 
does not address the anomaly problem; like related apporaches [2U [251 1261 12Z] , it rather 
avoids dealing with the problem by fixing the gauge or using deparameterization, choosing 
a time variable before quantization. With these additional steps, it is unlikely that the 
correct space-time picture is obtained in covariant form, and in fact the models evaluated 
so far have missed the signature change at high density. 

In this article, we introduce a new way of incorporating effects of inhomogeneity in loop 
quantum cosmology, dealing directly with wave functions. Adapting ideas of condensed- 
matter physics used to describe Bose-Einstein condensates]]] some effects of inhomogeneity 
will not be described by individual degrees of freedom but rather by non-linearity of wave 
equations for homogeneous models. The relationship between difference equations of loop 
quantum cosmology [291 US] and certain integrable non-linear Schrodinger equations has 
been noted in [31], providing additional motivation. The aim of the present article is to 
lay down the main ideas and to point out several new consequences for quantum cosmol- 
ogy. We find that tiny quantum corrections from inhomogeneous contributions to a large 
universe can add up coherently to produce sizeable effects on average, to be included in 
minisuperspace models. At this level, we will not yet address the anomaly problem and 
derive a detailed model, but we will demonstrate the prospects for this to be done at a 
later stage. 

2 Product states 

We start with a common way of dealing with inhomogeneity, viewing quantum space at 
a given time (a spatial slice used in canonical quantum gravity) as a collection of small 

1 Bose-Einstein condensates have been used in cosmological models before [55], but for matter rather 
than quantum geometry. Moreover, we are not suggesting that there is physical condensation, but rather 
use related mathematical techniques to describe approximately homogeneous geometries. 
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homogeneous parts. As one moves between spatial slices, the geometry evolves, resembling 
a many-body system of "interacting" elementary building blocks. Each building block 
(called a patch) has a quantum geometry described by a wave function of one of the 
well-known homogeneous models of quantum cosmology, and they all interact dynamically 
according to the quantized gravitational Hamiltonian. 



2.1 Classical model 

For simplicity, in this article we will assume isotropic patch geometries, determined classi- 
cally by a volume degree of freedom V^fc per patch, labelled by integers i, j, k to count 

patches in each spatial direction. A given spatial slice £ = \J^- k=l Vij t k — a differentiable 
manifold with a local atlas of coordinates — is then the union of M patches Vy,fc, or M 1 ^ 3 in 
each spatial direction. (Had we used anisotropic but still homogeneous patches, we would 
in general have three independent factors in M = A/1A/2A/3.) For now, we will assume 
M to be constant, which should be good for sufficiently brief evolution times. In more 
realistic models, the number M of patches should change in time, either by a fundamental 
process of discrete geometries being refined [321 [33], or by an approximation procedure 
akin to adaptive mesh refinement that maintains the decomposition into isotropic patches 
as a good model. (A time-dependent number of degrees of freedom is a general problem, 
studied for instance in [331 ESJ ESI E3 EE]-) 

For simplicity, we choose coordinates in space such that each patch has the same co- 
ordinate volume f v d 3 x = with & = Vq/M in terms of the total coordinate volume 
Vq of X (or of a large compact subset). The geometrical volume of each patch is then 
determined by the spatial metric which, if it is inhomogeneous, gives rise to different patch 
volumes Vij^- We assume that the metric is close to the one of a spatially flat, isotropic 
model with a longitudinal scalar mode, hob = a(t) 2 5 a b + 2L(t,x,y, z)5 a b- (We will use a 
lapse function corresponding to proper time, N = 1 — 2L/a 2 .) The patch volumes then 
take the values 

V iJtk = j dWdet h = a 3 / d 3 a;(l + 2L/a 2 ) 3/2 

f V 

3a / d 3 xL w — : + 3aL(x ijk )£ 3 (1) 



with the total volume V = a 3 Vo = ct'^N '. In the two approximations in the second line 
of this equation, we have first expanded the root and then replaced the patch-integrated 
L by its value at a point Xjj^ G Vi,j,&, such as the center. Since we assume the patches to 
be nearly isotropic and smaller than the variation scale of the perturbative inhomogeneity, 
L, both approximations are well justified. Solving (JTJ for L, we can therefore replace the 
continuum function L by deviations of the discrete variables Vij^ from the total volume 



V = Ma 3 t 3 : 



^).^. (2) 
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The dynamics of the Kj,jfc as functions of time is governed by a discretized version of 
the Hamiltonian constraint 

Hg TSV 4~ -^matter = (3) 

of general relativity, with contributions from the gravitational field and from matter. At 
this point, one will eventually have to face the problem of time and the anomaly problemU 
In this article, however, we focus on laying out the details of the new model, and therefore 
circumvent these difficult problems by formulating the dynamics in a specific gauge. With 
this choice, we may be blind to the complete quantum space-time structure, but new 
qualitative effects should still become visible. To proceed and to be specific, we assume 
matter to be dust, with Hamiltonian -ff ma tter = Pt/a 3 , where pt is a momentum variable 
conjugate to a matter degree of freedom t that will play the role of time. The role of time 
is made clear if we rewrite the Hamiltonian constraint equation as 

3 V 

Pt O Hgrav ~^Hg TSCV . (4) 

The variable p t then appears formally as an energy, or a canonical Hamiltonian that gener- 
ates evolution with respect to t. (More generally, we could assume matter to contribute to 
the Hamiltonian constraint by -ff ma tter — Pt' / a 3( - 1+M ^ if there is a perfect fluid with equation- 
of-state parameter w. A time variable t' different from t then parameterizes evolution.) 

To derive the dynamics in detail, we start with the classical Hamiltonian constraint 
of general relativity and write it in discrete canonical variables together with their 
momenta n^fc, related to Vij^- 111 the ADM formulation of canonical gravity, the spatial 
metric h ab is canonically conjugate to 

^ = ^*h {K a b _ K c h a b) (5) 

lo7rG 

(See jl5] for an introduction to canonical gravity.) We compute the canonical variables 
in our perturbed situation by writing h ab = h5 ab + 5h ab and n ab = (7i/V )S ab + 5ir ab , split 
into background variables h = a 2 and 7? (spatial constants) and inhomogeneity 5h ab and 
57i ab . We divide 7? by V in n ab to ensure that the symplectic term J s d 3 xh a i ) iT ab = tin + • • • 
assumes the canonical form in its background term. For scalar modes in longitudinal gauge, 
5h ab = 2L5 ab and 5n ab = 5n5 ab . 

To avoid overcounting of degrees of freedom, we require the inhomogeneity 5f of any 
field / = / + Sf to satisfy L d 3 a;5/ = when integrated over all of space. (We turn 
inhomogeneities of tensor fields such as 5h ab into scalars using the background metric 8 ab .) 
As a consequence, / = J s d 3 xf is indeed the spatial average. At this stage we do not assume 
that 5f is of first or any specific order in perturbation theory; we have simply rearranged 
our degrees of freedom by splitting them into background variables and inhomogeneity. 



2 Regimes of interest here, with small inhomogeneity, are usually semiclassical regarding quantum ge- 
ometry. In such a situation, effective constraints [39l [40l [41] can be used to solve the problems of time 
021 gll |44] and anomalies [3]. 
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The symplectic structure, our current interst, only refers to degrees of freedom but not to 
orders of perturbation theory: higher perturbative orders do not introduce new degrees of 
freedom. We will introduce the perturbative expansion when we prepare our Hamiltonian 
for a derivation and analysis of equations of motion. 

Any terms linear in Sh a b or 5v: ah in the Hamiltonian or symplectic term J s d 3 xh a b7i ab 
therefore vanish. Inserting the inhomogeneous metric h a b = (a 2 + 2L)S a b in (jSJ) (with 
vanishing shift and longitudinal lapse in K a b) results in the inhomogeneous momentum 

* ab = -^(* + *(^Y)t ab (6) 



8nG V W, 

from which we read off the momentum of h = a 2 as tx = — dVo/87rG with the total coordi- 
nate volume Vo, and the momentum of 5h a b = 2L5 a b as 5ir ab = — (a(L/a 2 )' /8iiG)5 ab . By a 
canonical transformation we can switch to volume variables as defined in our patch model: 
we have momenta 



Tl v = and IL k = — ^ (7) 

V 12nGV hhk 12nG \ V J K J 

of V and Kj,fe- 

For small inhomogeneity, it is sufficient to expand the Hamiltonian constraint to second 
order in 5h a b (or L) and its time and space derivatives. Starting from 

# g rav = J d " xN (K^Kab ~ K 2 - 3 R) Vdrth = J d 3 xH gTav , (8) 

we obtain 

3 ( , 2 L 2 — A(g/a)LL + 4(d/ d) 2 L 2 _*({ dL\ 2 , , u » 
?W~-^^a + - +a ^ \^— b j - -L ^ 

In the Hamiltonian -ff grav = J s d 3 x"H grav we can integrate by parts in spatial derivatives, 
replacing second-order derivatives by first-order ones. (Boundary terms will play no role 
in what follows.) Moreover, it turns out that the time derivatives of L can be written more 
compactly if we use L/a 2 , a combination of variables that is also more convenient when 
expressed by patch volumes: L 2 — 4(a/a)LL + 4(a/a) 2 L 2 = ((L/a 2 )') 2 . The Hamiltonian 
density we use will therefore be 





*W = szf - + -2 +^2J-So-^ • ( 9 ) 



We then introduce our background momentum a/a — — AirGTly and the patch momenta 
(L/a 2 )' — AirGHij^ after replacing the integral by a sum over patches, J s d 3 x?/ grav 
Si j k%j>k- O ur discretized Hamiltonian then is 



= ~^ GV fa + jfH K,k + ■■) (io) 
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where the dots indicate the derivative terms after discretization. 

We have quadratic single-patch Hamiltonians in the first two terms, analogous to har- 
monic 1-particle Hamiltonians of our many-body problem. Spatial derivatives of L must be 
discretized before they can be expressed in terms of the Vij^- The discretization procedure 
is a matter of choice and, to some degree, convenience; we will make use of 



Quadratic expressions of spatial derivatives in ([9]) then provide interaction terms that 
can be written as depending on either the patch geometries in product form, such as 
V (i,j,k)+b V (i,j,k)-h or more conveniently, the difference (V {i>j>k)+ ~ b - V^j.j) in discrete min- 
isuperspace. The latter version is closer to interactions of many-body systems depending 
on the distance between particles. 

In addition to interactions between neighboring patches, each patch volume interacts 
with the average volume V because it appears in some factors in the Hamiltonian. These 
variables are not independent but satisfy Yliijk^hi^ — V. In order to focus on the self- 
interaction of inhomogeneity, we will treat V as an external parameter for the dynamics 
of the Vij t k, corresponding to the common approximation in cosmology that ignores back- 
reaction of inhomogeneity on the background. 

2.2 Quantization 

Each patch of volume V^j^ and expansion rate related to Hij,k is isotropic and may be 
quantized as a single minisuperspace model, corresponding to the 1-particle Hilbert space 
of a many-body system. One may follow either Wheeler-DeWitt quantization or loop 
quantization, both with volume representations in which becomes a multiplication 
operator. In the former case, one deals with wave functions ^(^ij.fc) m L 2 (WL + ,dVij t k) 
and the momenta act by flij,k = —ihd/dVi^^; in the latter, ipv ijk is an element of the 
non-separable sequence space £ 2 (M.) and exponentials of momenta, rather than momenta 
themselves, are quantized: 



action of exp^tfjj^ITjj^/ft) on the sequence space is not continuous in 5ij,k, and a derivative 

3 Choices in the discretization procedure affect physical implications and can therefore be tested for 
their consistency. These choices are related to quantization and regularization ambiguities in canonical 
quantum gravity from which our expressions should follow in some complicated way. If discretization 
choices can be restricted, the same will be true for ambiguities of an underlying fundamental theory. 

4 We work in a special Abelian sector of homogeneous loop quantum cosmology. The general non- 
Abelian structure is more complicated due to refinement features [35], but qualitative aspects are shown 
well by the Abelian simplification. 





indicating by b the unit vector in the 6-directionJf] 
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by Si,j,k, which would otherwise result in an operator for ITj^, does not exist. (A second 
difference between the quantizations is that Kj.fc i n Wheeler-DeWitt models is usually 
taken as the (positive) volume, while loop quantum cosmology is based on triad variables 
in which Vij } k is the oriented volume, which can turn negative if the orientation is reversed. 
We therefore use the full real line R in the sequence space, rather than R + . Note that this 
resolves self-adjointness issues of derivative operators on L 2 (M + , dVij^)-) 

Both representations are well-defined but not unitarily related to each other; they 
lead to different physics. Especially at high curvature, where H-ijk is large, effects of the 
loop quantization can differ significantly from those of the Wheeler-DeWitt quantization. 
The discreteness inherent in shift operators ( 1T2|) relating derivatives is then important, in 
addition to the discreteness implemented by our treatment of inhomogeneity. 

If inhomogeneity is small, the patches evolve nearly independently of one another with- 
out strong correlations, and the evolved state remains a product state \&(Vi, V 2 , . . .) = 
^1(^1)^2(^2) • • • of the individual patch wave functions ipi if the initial state is of such a 
form. Each single-patch wave function evolves according to a differential (Wheeler-DeWitt 
[17] ) or difference (loop quantum cosmology [29l 150] ) equation if inhomogeneity can be 
ignored. With inhomogeneity included, interaction terms between the individual wave 
functions occur on superspace, complicating the dynamics. If inhomogeneity is sufficiently 
small, however, the interactions can be treated by approximation, such as perturbation 
theory. 

Small inhomogeneity at the level of quantum geometry also implies that the individual 
wave functions are very similar to one another, so that the full state can approximately be 
written as \P(Vi, V2, ■ ■ •) = ^(^1)^(^2) ' ' ' with a single wave function ip to be solved for. 
This form of product states allows one to map many-body dynamics to 1-particle dynamics 
in a specific potential, described by a wave equation that turns out to be non-linear. At 
this stage, standard techniques to describe matter condensates, in which individual wave 
functions of different particles are exactly equal to one another, can be applied. 



2.3 Condensate 

By our preceding considerations in cosmology, we have realized a mathematical formulation 
with all the ingredients used in the description of Bose-Einstein condensation. We interrupt 
our discussion of cosmology to recall salient features of this important system in condensed- 
matter physics. In this example, \l/ is a many-body state, and ip the 1-particle wave function 
common to all constituents of the condensate. Taking the same ip is n °t an assumption 
because condensed particles have exactly the same wave function. 

Assuming pointlike interactions between the particles, described by a delta-function 
potential of strength a, we have the many-body Hamiltonian 

k = E f + v w) + \ a E - *i) ( 13 ) 

i=i v / #j 

for n particles of mass m in individual potentials V(xi). With a product state ^(xi, x 2 , ■ ■ •) = 
ip{x\)ip{x2) ■ ■ ■ for the condensate, we compute the expectation value of the Hamiltonian 
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as 




(14) 



The first term just adds up the 1-particle expectation values computed for the wave function 
if). The second term is not equal to a 1-particle expectation value. However, we can 
formally interpret it as the expectation value of a "potential" |?/>(x)| 2 depending on the 
wave function. Accordingly, the 1-particle dynamics and energy spectra are governed by a 
non-linear Schrodinger equation, the Gross-Pitaevski equation 



For a full and rigorous derivation, see [1B1 US]- 

Interacting many-body dynamics of the condensate wave function can therefore be 
mapped to non-linear 1-particle dynamics. 

3 Non-linear dynamics in quantum cosmology 

With the preparations presented in the preceding section, we propose a new method to 
deal with small cosmological inhomogeneity, making use of the same ideas and initial 
constructions employed to describe matter condensates. Well-established methods then 
provide a tractable approximate description by non-linear dynamics of a homogeneous 
model. 

3.1 Equation of motion 

Except for the differences in the conceptual nature, regarding for instance the approxi- 
mations and assumptions used, our model for inhomogeneous quantum cosmology so far 
resembles those of matter condensates rather closely. The main mathematical difference 
lies in the interaction potential. For particles in a condensate, a delta function of the 
distance between particles is a good approximation for nearly pointlike interactions, which 
can be smeared out to more-complicated functions for realistic systems. The interaction 
potential we obtain in cosmology, expanding and discretizing the gravitational Hamilto- 
nian constraint, is a quadratic polynomial in the distances in minisuperspace. Although 
the single-patch wave equation we obtain is still non-linear, as in the presence of any kind 
of interactions, it is more complicated than in the Gross-Pitaevski equation. 

Another difference between the models is the discreteness of the quantum representation 
used in a loop quantization, in addition to the discretization of space by patches Vjj,*.. 
Not only space but also superspace is then discrete. As a consequence, wave equations 
in loop quantum cosmology are difference equations, and with our method to include 
inhomogeneity we will be dealing with some version of a discrete non-linear Schrodinger 
equation, one example given by 




(15) 




(16) 
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However, since we are not dealing with pointlike interactions in superspace, modeled by 
delta functions, but rather with polynomials, the non-linearity will be different. In fact, 
our equation will not only be non-linear but also non-local but nevertheless, as it turns 
out, well-suited to canonical effective methods. 

Using the same starting point as in Bose-Einstein condensation, the key step is to 
evaluate the expectation value of the interaction Hamiltonian in a product state. To 
illustrate the main consequence, we consider just two variables V\ and V 2 interacting with 
each other via a potential Wi n t(Vi, V 2 ) = ce(Vi — V" 2 ) 2 /V 2 as in a discretized (Q. We divide 
by the total volume squared, treated as an external but time-dependent parameter, in 
order to have the correct scaling behavior of the Hamiltonian under a change of the spatial 
region. The expectation value of the quantized W int then produces a term 



where we introduce 8V := V2 — V\. 

We can perform the second integration independently of the first over V\. It depends 
on the wave function, but if we assume that ip is sharply peaked around the expectation 



equals the quantum fluctuation of V in the state t/>(V). Instead of a non-linearity potential 
depending on ip{V) or ip n as in ( fl6|) . we have a non-linearity potential that depends on the 
wave function via moments such as AV^. For instance, following the preceding arguments 
and noting that the minisuperspace V is quantized to a discrete parameter n, we need to 
consider an equation of the form 



We note that equation (1201) is not only non-linear but also non-local: the coefficient 
(An) 2 = ^2 n (n — {n)^) 2 \ijj n \ 2 depends on all values of ip n . Moreover, the equation as 
written is meaningful only for n 7^ 0. At n — 0, the volume vanishes and we encounter 
a cosmological singularity. By inverse-triad corrections [50], loop quantum cosmology re- 
solves this singularity in such a way that 1/n is replaced by a bounded function. For 
simplicity, we will not discuss these terms here and instead focus on evolution at large n. 

We must ensure that our assumption of a sharply peaked state remains true for the 
approximation to be valid. If the state is not sharply peaked or if the approximation is to 
be driven to higher orders, we can use a derivative expansion of if). Writing 




(17) 






(20) 



mV 1 + SV)\ 2 = \ij((V)+5V + (V 1 



(V)))\ 2 
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and expanding by V\ — (V), we obtain 



/ 



dV 1 \ip{V 1 )\ 2 W oa na a (V 1 ) 



with the non-linearity potential 



^ non i in (y) 



T,^vy pU) (v-(v}) 



3=0 



(21) 



where the moment (AV") 2 ^) is the V^-fluctuation computed with the "distribution" pv\ 
defined as the j-th derivative of p(V) = \i[)(V)\ 2 . Note that these derivatives need not be 
normalized or positive, so that we do not have probability distributions and fluctuations 
in the statistical sense. Nevertheless, the resulting numbers are well-defined as parameter- 
izations of the non-linearity potential. 

We continue with a discussion of the leading-order equation fl20l) . 

3.2 Solution procedures 

An inverse scattering transform is the method of choice to solve non-linear discrete or 
differential Schrodinger equations [51]. However, the equation we obtain here, (!20l) . is not 
only non-linear but also non-local. Standard techniques are therefore not readily available. 

Non-local equations can sometimes be treated by replacing the non-local coefficient by 
new auxiliary degrees of freedom, as in [52] in the context of the non-linear Schrodinger 
equation. If the new degree of freedom is subject to a differential or difference equation 
with a source term given by the original wave function t/>, its general solution is a non- 
local expression in ip (integrating its product with the Green's function of the auxiliary 
equation). If the right equation is chosen, the general solution for the auxiliary variable 
may provide the non-local coefficient, (An) 2 in our case. Here, however, such a treatment 
is not obvious. 

Instead, canonical effective methods [53] based on the dynamics of moments of a state 
provide solution techniques well-suited for equations such as (I20I) . The non-local coefficient 
is a second-order moment of the wave function; using equations for the moments instead of 
■0n itself then provides a reformulation of the problem in variables in which the non-locality 
disappears. Morally, this procedure is a version of introducing new degrees of freedom 
related to the wave function non-locally, for moments^ such as A(n a ) := Y2 n ( n ~ ( n )) a |'0n| 2 
with the expectation value (n) = ^2 n n\ip n \ 2 are non-local in ip n . However, in quantum 
physics the moments are not auxiliary variables but rather variables of prime physical 
interest. For a = 2, we have quantum fluctuations, and higher moments with a > 2 
provide additional statistical information about the state. 

For linear discrete or differential Schrodinger equations, canonical effective techniques 
[53] amount to a systematic expansion of Ehrenfest's equations, used not just to derive the 

5 Our notation is a variation of the common (An) 2 = A(n 2 ) at second order, the parenthesis displaced 
to be unambiguous at higher orders. 
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semiclassical limit in rigorous terms [51] but also to compute quantum corrections to any 
desired order in h. For our purposes, we need to generalize these methods to non-linear 
equations as encountered here. 

In quantum mechanics, a set of iV basic operators Jj with closed linear commutators 

[J, Jj] = '^^Cij k Jk (22) 

k 

(perhaps including the identity operator if some commutators are constants) provides a 
closed algebra for expectation values under Poisson brackets 

{<4>,<4>} = %^- (23) 

If the operators are complete, any observable can be expressed as a function of the expec- 
tation values ( Jj) and moments 

\ i / \ i I symm 

with operator products in totally symmetric ordering. Using linearity and the Leibniz 
rule for Poisson brackets, these expectation values and moments form a Poisson manifold. 
Their dynamics is determined by the Hamiltonian flow generated by the expectation value 
Hq := (H) of the Hamiltonian constraint, another observable interpreted as a function of 
expectation values and moments. Hamiltonian equations of motion usually couple infinitely 
many moments to the expectation values, but a semiclassical expansion to some finite 
order in fi results in finitely coupled equations which can be solved at least numerically. 
Computer-algebra codes exist to automate the generation of equations to rather high orders 
[55] (so far restricted to canonical commutators). 

Writing Jj = (Jj) + (Jj — (Jj)) in the quantum Hamiltonian Hq = (7J(Jj)) and per- 
forming a formal expansion in (J — (Jj)), the Hamiltonian flow is generated by 

„ Q = „ m) ^-...-— —j—^Jl^ j . (25 , 

The first term is the classical Hamiltonian evaluated in expectation values, and the series 
includes quantum corrections of progressing order £\ aj. Equations of motion follow from 
Poisson brackets. 

These constructions rely on commutators of linear operators and cannot be used directly 
for non-linear Schrodinger-type equations. Nevertheless, a closely related procedure can 
be followed for equations such as (120]) in which the non-linearity comes from non-local 
coefficients depending on the moments. As one can readily confirm by computing time 
derivatives of expectation values directly using (12"U|) for wave-function factors, the evolution 
of moments is now governed by a quantum Hamiltonian (|25l) in which one initially treats 
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the moments that appear in the non-local coefficients as external functions; because they do 
not come from a linear operator, they do not appear in commutators or in Poisson brackets 
of the moments when equations of motion are derived. In the equations of motion, once 
derived, these variables are to be equated to the moments they signify, providing additional 
coupling terms between moments compared with a linear Hamiltonian. 

In our case, the Jj are given by three basic operators, a multiplication operator by n 
(the volume operator) and two shift operators h and h' that change n by ±1, implementing 
ffl2|) with 8—1. In terms of canonical variables (n, P), we can write shift operators as 
quantizations of h = exp(iP). The commutators 

[n, h) = -hh , [h, h ] ] = hh ] , [h, h ] ] = (26) 



then define the basic algebra ([22]) of our loop-quantized theory, and correspondingly the 
Poisson brackets of expectation values and moments of n and h. Moreover, since we 
introduced complex variables, the reality condition hh = 1 as well as analogs for the 
moments (such as A(hh*) = 1 — hh*) must be satisfied. 

We can realize the linear part of (|20|) as the Schrodinger equation with Hamiltonian 
operator h + fv — 2, resulting in the quantum Hamiltonian Hq 1 = h + h* — 2 depending only 
on expectation values but not on moments. Adding the non-linearity, we have an extra 
term — ^aA(n~ 2 ) with A treated as a constant to be set equal to A = (An) 2 in equations 

of motion, and (n~ 2 ) to be expanded by moments as in fl25|) . For the difference equation 
(120]) . we then have the quantum Hamiltonian 

H Q = h + h*-(2- aAip 2 )) = h + h* -2 + aA{3n~\An) 2 - 20n~ 5 A(n 3 ) + • • •) (27) 

with the non-local coefficient A treated for now as an external parameter. (Instead of the 
inverse of n, which is ill-defined at n — 0, modifications due to inverse-triad corrections in 
loop quantum cosmology should be used at small n [50].) 

We obtain the equations of motion from Poisson brackets, in which we then set A = 
(An) 2 : 

h = i(h-h*) (28) 

h [AnV (An) 2 A(nh) . . 

h = \2ia- J - 6ia- — - + ■■■ (29) 

n \ n J n 4 

2 

= 2i{A{nh) - A{nh*)) (30) 



d(An) 



dt 

and so on for further moments. Instead of n = as in the linear case, we can combine the 
first two equations to obtain 

• • h + h* /An\ 4 (An) 2 (A(nh) + A(nh*)) , , 

n = i(h - h*) = -12a— — — - 6a± — J -— — ^ + ■ ■ ■ . (31) 

n \ n J n 4 

Non-zero moments imply (negative) acceleration of the volume expansion. 
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3.3 Interpretation 

Irrespective of the precise form of non-linearity, its presence has several general conse- 
quences of potential importance for quantum cosmology. An obvious and seemingly prob- 
lematic implication is a loss of unitarity: wave functions evolved by the non-linear equation 
do not have preserved scalar products with other evolved states. There is no linear op- 
erator that could serve as a Hamiltonian whose adjointness properties one could analyze 
by standard techniques. Still, a straightforward direct calculation shows that the norms 
of states (but not scalar products (4>\ip) of different states) are preserved. How- 
ever, the original many-body system is clearly unitary, and therefore non-unitarity is a 
consequence of the reductions and approximations used. In order to interpret the non- 
linearity correctly, we should therefore look back on the constructions used to descend 
from many-body dynamics to a 1-particle equation. 

For a matter condensate, we obtain the non-linear wave equation ( fl5l) in a rather 
indirect way: We do not reduce the many-body wave equation for \l/ directly, but rather 
compute the expectation value of the Hamiltonian (TH1) . rewrite it in terms of the 1-particle 
wave function if>, and recognize the extra term as a formal analog of a potential depending 
on the wave function. This potential, inserted in the standard Schrodinger equation, then 
provides ffl5|) . a step which is again only formal. Experience shows that the resulting non- 
linear equation nevertheless captures crucial properties of the many-body problem, and 
rigorous proofs have been provided [HI EH] . 

One can avoid the last formal step by forgoing wave equations and instead using the 
expectation value ( !T4|) to compute the spectrum of the many-body Hamiltonian, for in- 
stance by variational methods applied to the 1-particle wave function if) on the right-hand 
side of ( fl^l) . If the spectrum of the Hamiltonian is known, evolution properties then follow 
without directly using the non-linear equation ( fl5l) . Similarly, effective canonical equations 
in quantum mechanics refer to expectation values of the Hamiltonian, such as (fT4|) rather 
than wave equations, and are therefore less sensitive to the apparent loss of unitarity. 

The physics of the system therefore does not suffer from a lack of unitarity. Moreover, 
since the norm is still preserved, the probability interpretation of a single state remains 
meaningful. Instead of using f lTo^) as a fundamental wave equation for some function if> in a 
Hilbert space, the equation models other dynamical effects, such as the evolution of particle 
distributions or the approach and possible interaction of superposed states. Properties 
such as the overlap of superposed states or the distance between different distributions 
can be determined from moments of a single wave function for the superposition and 
are independent of scalar products of the wave function with other states; they can be 
analyzed with a formal equation lacking unitarity. These are also the properties that 
effective equations are sensitive to. In quantum cosmology, such questions are usually of 
most interest because the exact state or wave function of quantum space is not accessible 
by observations available now or in the foreseeable future. Our model and with analog fj20|) 
of the Gross-Pitaevski equation (fTS]) is therefore reasonable. 
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4 Discussion 



We have introduced a new model for inhomogeneous quantum cosmology, aiming to cap- 
ture essential features of the interacting dynamics of different parts of quantum space. The 
processes we describe therefore provide the dynamics of structure formation at a funda- 
mental level. Using several approximations, justified when inhomogeneity is sufficiently 
small, and importing ideas of condensed-matter physics, we have been able to map the 
complicated many-body dynamics to a non-linear minisuperspace equation^] 

In addition to the approximate nature, several differences with the condensate model 
occur: 

• In the cosmological model, "interactions" between different patches are realized in 
superspace, not in actual space. Patches do not interact depending on their spatial 
distance, but depending on what their geometries are: The gravitational Hamiltonian 
depends on inhomogeneous modes, or on deviations of patch geometries from the 
spatial average. 

• There is no delta-function potential (for pointlike interactions) but rather a polyno- 
mial potential, obtained by expanding the gravitational Hamiltonian as a function of 
patch geometries. As a consequence, the non-linearity is realized non-locally in the 
configuration space of wave functions. 

• While the many-body Hamiltonian of a condensate is well known but difficult to 
deal with, a consistent version of an inhomogeneous gravitational Hamiltonian in 
quantum gravity is still lacking. In particular, covariance conditions and the related 
problem of anomalies have not been evaluated in sufficient detail [601 ED E21 E21 EI] . 
(But see [651 ESI EH ES EH] for recent progress.) 

In this situation, having an approximate description of incompletely known dynamics, 
we cannot expect to derive detailed quantitative cosmological scenarios. (This statement 
does not only apply to our new method, but to all derivations possible in quantum cos- 
mology so far.) Effective techniques, as used in our solution procedure for non-linear 
non-local equations, provide means to parameterize ambiguities and ignorance, and to dis- 
cuss anomalies, but no details are available yet. We therefore focus our discussion on new 
qualitative features suggested by the non-linearity of the homogeneous model. 

Non-linear wave equations provide new forms of minisuperspace effects that capture 
crucial properties of averaged inhomogeneity. These terms need not require high, near- 
Planckian densities to be significant because they could potentially be large when many 
patch contributions are added up, even if each of them is tiny. All leading contributions 
have the same sign because they come from volume fluctuations, required to be positive. 
No cancellations happen when one sums over all patches, potentially giving large effects. 

6 Other versions of non-linear quantum cosmology have been proposed [56j EH ES [59], motivated by 
non-commutativity and information-theoretic arguments. 
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For certain behaviors of quantum fluctuations as functions of time or the volume, our non- 
linearity can be interpreted as a cosmological-constant term, which turns out to be negative. 
(Again, the sign is determined because quantum fluctuations are always positive.) It 
remains to be seen whether more-refined models, including those with anisotropic patches, 
or higher orders in the moments in (12ip . not all of which are restricted by positivity, as well 
as perturbed Hamiltonians beyond second order can turn the sign to provide an overall 
positive cosmo logical constant. 

An interesting feature of non-linear wave equations is the existence of a particular type 
of solutions: solitons. These are sharply peaked wave packets which evolve without chang- 
ing shape. Moreover, if solitons occur in superposition, moving in different directions, 
they may occasionally overlap but do not influence each other. After they have moved 
through the same spot, they retain their old shapes. Such states are a promising candi- 
date for new dynamical coherent states in quantum cosmology. In contrast to kinematical 
coherent states (or Gaussians) commonly used in such cases, solitons are adapted to the 
dynamics and, in the indirect way that employs non-linear wave equations, capture prop- 
erties of inhomogeneity; in fact, their existence as solutions relies on deviations from exact 
homogeneity. 

The existence of solitons and the integrability of equations, together with the associ- 
ated possibility of chaos, depends sensitively on the form of discrete equation [70]. The 
discreteness, in turn, is related to quantization and regularization ambiguities in canonical 
quantum gravity. The strong sensitivity of some physical features may allow one to find 
tight restrictions on ambiguities. 

We end by mentioning another, more speculative consequence. In quantum cosmology, 
solitons in superposition would correspond to different universes superposed in the same 
state. Solitons may overlap but do not affect each other's motion; they always form 
separate contributions to the total state. Solitons and the non-linear wave equations they 
solve could therefore play a role in the description and analysis of multiverse models. 
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